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Effects of magnetic field on the quadrupolar ordering are investigated with inclusion of fluc- 
tuation of order parameters. For the simplest model with the nearest-neighbor quadrupolar 
interaction, the transition temperature and the specific heat are derived by the use of the re- 
cently proposed effective medium theory. It is shown that magnetic field H has two competing 
effects on the quadrupolar ordering; one is to encourage the ordering by suppressing the fluctua- 
tion among different components of order parameters, and the other is to block the ordering as in 
antiferromagnets. The former is found to be of order H 2 and the latter of order H 4 . Hence the 
fluctuation is suppressed for weak fields, and the transition temperature increases with magnetic 
field. The fluctuation effect is so strong that the entropy released at the quadrupolar ordering 
is only about half of the full value In 4 even without the Kondo effect. 
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§1. Introduction 

Some compounds with orbital degeneracy show order- 
ing of quadrupole moments without magnetic order. A 
typical example is CeB6 which also shows the Kondo ef- 
fect with the Kondo temperature T K ~ 1KB The Ce 3+ 
multiplet (4f\ J = 5/2, S = 1/2, L = 3) is split un- 
der the cubic crystal field into the ground-state quartet 
Tg and an excited-state doublet T-? separated by about 
540K.EP It has a transition from the phase I (paramag- 
netic phase) to the flkase II at 3.3K under zero field. 
NeutronB and NMRfJtP experiments show that electric 
quadrupole moments are ordered with a periodic pat- 
tern in the phase II. Hence it is also called the antiferro- 
quadrupolar phase (AFQ). With further decrease of tem- 
perature, another transition occurs at 2.4K to the phase 
III, which accompanies the magnetic ordering with a 
complicated periodic pattern. The phase boundary be- 
tween the phases I and II shows unusual dependence 
on magnetic field; the transition temperature Tq(H) in- 
creases with magnetic field H. Namely, while a mag- 
netic field will ordinarily destroy an ordered phase, it 
assists the quadrupolar ordering in CeBg. The phase 
boundary shifts tnJower temperatures in the dilute al- 
loy Cei.zLazBgBEP The AFQ is found also in TmTeJf 
which is insulating. Though its phase diagram has some- 
thing similar to that of CeBg , Tq (H ) shows a maximum 
as a function of magnetic field in contrast to that of 
CeBg. In the latter, Tq(H) continues to increase up to 
the maximum field measurable at present. 

One of the keys to understanding the transition be- 
tween the phases I and, II should be the result of specific- 
heat measurement pt3 which indicates large fluctuation 
there. Although the peak at the AFQ transition is small 
in zero magnetic field, it grows larger with increasing 
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magnetic field. The small peak suggests that only a small 
change of entropy is involved at the transition, and that 
a large amount of short-range order should exist even 
above Tq. In other words, the long-range order is dis- 
turbed by some large fluctuation. Therefore, the experi- 
mental result of specific heat suggests that the large fluc- 
tuation should be suppressed by a magnetic field. The 
purpose of this paper is to clarify how a magnetic field 
affects the fluctuation at the quadrupolar ordering. We 
take account of fluctuation in our calculation with use of 
an effective medium theory for quantum spins and mul- 
tipoles.tll 1 

Several models for the quadrupolar ordering in CeBg 
have been proposed. Hanzawa and KasuyaliJ assumed 
that the ground state is the doublet Tt, and consid- 
ered the quadrupolar interaction among Ce ions. In 
this model, large off-diagonal elements of quadrupole- 
moment operators between IV and T$ become impor- 
tant as magnetic field increases. Thus, Tq(H) increases 
with H in the mean-field approximation. Although this 
model should be effective for the system with the T-? 
ground state, it was shown later by experiment that Tg 
is the ground state and T7 is an_.expitcd staterP An- 
other model proposed by Ohkawali3't3 has the symmet- 
ric form between the spin doublet and the orbital dou- 
blet represented by pseudo spins. The Hamiltonian in- 
cludes an intersite interaction which is biquadratic in 
spin and quadrupole moments. This interaction en- 
hances the quadrupolar interaction with increasing H. 
Hence Tq(H) increases in the mean- field approxima- 
tion also in this model. Shiina et.alS3 recently rewrote 
Ohkawa's model in the form including magnetic oc- 
tupoles and dipoles, and solved the resultant model by 
the mean-field approximation. Finally, Uimin ancLl, 
present authors treated the effect of fluctuation.!!^ I 
They approximated the pseudo-spins by classical spins 
and applied the spherical model which can include flue- 
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tuation in a simple manner. As a result, the transition 
temperature increased with magnetic field in contrast to 
the result of the mean-field approximation. 

In this paper we accomplish substantial improvement 
over the theory of refjl^. The fluctuation emphasized 
there is due to the competition among order parameters 
with different wave numbers k. The intersite interaction 
was taken to be the electrostatic one between the IV 
type quadrupole moments. Strength of the quadrupolar 
interaction at k = (tttttt) is almost the same as that at 
k = (irirO) in the case of the electrostatic interaction. 
Competing fluctuations are suppressed by a magnetic 
field. However, as will be discussed in detail later, what 
is essential of the magnetic field is to suppress compe- 
tition among components of fluctuating order parame- 
ters, namely, O 2 and 2 2 . The presence of competi- 
tion among order parameters with different wave num- 
bers helps the suppression effect mentioned above. The 
value of dTq (H) / dH for the electrostatic interaction is 
about twice as large as that for the nearest-neighbor one. 
Although this enhancement is substantial quantitatively, 
we regard that the specific form of interaction is not es- 
sential for qualitative understanding. Hence we adopt 
the nearest neighbor interaction for simplicity of calcu- 
lation in this paper. Further, we do not restrict the 
quadrupolar interaction to that of the I^-type moments. 

§2. Model for Quadrupolar Interaction 

2.1 Model 

Our model corresponds to a localized electron system 
which has interaction between quadrupole moments only. 
Possible targets of the model are not only to CeB6 but 
also related quadrupolar systems under the conditions 
to be described below. An insulating compound TmTe 
has a phase diagram similar to CeBg. If the increase 
of transition temperature in these two materials comes 
from a common mechanism, the interaction with conduc- 
tion electrons should not be essential. Hence we do not 
consider this interaction explicitly. 

Cerium sites in CeBg constitute the simple cubic lat- 
tice. At temperatures near Tq, the excited C-EF states, 
which are separated by about 540K in CeB6,EP have lit- 
tle influence. Therefore we deal with the ground-state 
quartet Tg only. We represent here the basis set of the 
T$ states with use of the \J Z ) basis, 



i± 



±- 



2±; 



1 



±- 



(2-1) 

The Hamiltonian for our model is composed of the 
quadrupolar interaction term Wj n t and the Zeeman term 
Hz: 

H = H int + J2 n z(i)> (2-2) 

i 

H int = Jrs E [O2 (i)O 2 °(j) + 2 2 {i)0 2 2 {j)} 

(ij) 

+J r5 ^ [Oy Z (i)O yz {j) + O zx (i)O zx (j) + O xy (i)O xy (j)} , 

(ij) 

(2.3) 



(2.4) 



where 6/7 in the Zeeman term is the g factor. The sum- 
mation in 7ii n t is performed over all the nearest neigh- 
bor pairs. Quadrupole-moment operators at each site 



are given as 
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2 "= [3(J^- J(J+1)], 
O2 = {J x f - {J y )\ 



o yz = j y j z + j z j y , 



o zx 

X v 



J Z J X 
j x j y 



j x j z , 
j v j x . 



(2.5) 



(2.6) 



The interaction constant Jr3 is different from Jrs in 
general. We take three types of interaction constants: 
(i) T3-type interaction Jr3 > 0, Jrs = 0; (ii) Ts-type 
interaction Jr3 = 0, Jrs > 0; (hi) (r3 + rs)-type inter- 
action 16Jr3 = Jrs > 0. In the case (hi), interaction 
constants are taken so that every ordering of five types 
of quadrupole moments has the same transition temper- 
ature in zero magnetic field. This becomes clear in the 
pseudo-spin representation described in the next subsec- 
tion. It is expected that the effect of fluctuation is the 
largest in the case (hi). 

This model is the simplest one for the quadrupolar or- 
dering in CeB@. However, in the literature the increase 
of the transition temperature has never been discussed 
within this model. This is probably because the result 
of the mean-field approximation does not show the in- 
creasing transition temperature. We treat this model 
beyond the mean-field approximation to discuss effects 
of the fluctuation. 

2.2 Pseudo-spin representation 

It is convenient for calculation of Tq(H) to use the 
pseudo-spin representation. When the Hilbert space is 
limited to the Tg states, all operators are represented by 
these pseudo-spins because all the matrix elements can 
be represented by those of pseudo-spins. The pseudo- 
spin operators are defined by 

t z \^i±) = \A±), ry 2± ) =-|Va±>> 

a z \ij a+ ) = \^ a + ), £T*|Va-) = -|Va->» (2-7) 

a + \ip a -) = \i> a +), o--\i> a+ ) = \lp a -), 

(a = 1,2). 

Note that this definition is different from that in refs.|l2], 
|l4] and |l^ where the positive eigenvalue of the pseudo- 
spins is set to 1/2. It is set to unity here to remove the 
constant 1/2 from various expressions. 

In this representation, quadrupole-moment operators 
are given by 



n ® ^ 2 n 2 ^ ^ 



V3 



X y — 



V3 



(2.8) 



(2.9) 
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The interaction Hamiltonian is rewritten as 

Hmt - Jf 3 E b-'W^O") +r x (i)r x (j)} 

(ij) 

+JUY, Ty WT y (j)<T(i)->*U), (2-10) 

(ij) 

with Jp 3 = J r3 x 64/3 and Jf 5 = Jrs x 4/3. Hence 
Jf 3 = Jf 5 is satisfied for the (T 3 + T 5 )-type interaction. 
The Zeeman term is rewritten as 

Wz(») = - E ^{l + ^T a {i)\^{i)H-, (2.11) 

a=x,y,z ^ ' 

where T Q = -±={3( J Q ) 2 - J (J + 1)} x ^ i s g i ven i n the 
pseudo-spin representation by 

T X = --T Z + —T X , T y =--T Z -—T X , T Z =T Z . 

2 2 2 2 

(2.12) 

Intuitively, (1 + 4/7 T") represents the magnitude |m| 
of the magnetic moment with 3/7 < |m| < 11/7 and a a 
represents the orientation of it. 

§3. Static Approximation for the Effective 
Medium 

We recently formulated a dynamical effective medium 
theory for quantum spins and multipoles.EIP The theory 
is a generalization of the spherical model approximation 
for the Ising model.li§liif It takes account of fluctuation 
up to 0(X/z n ), where z n is the number of interacting 
neighbors. We use the static approximation (SA) which 
neglects the dynamical nature of fluctuations. Although 
it is difficult to estimate the accuracy of the SA in the 
temperature range of the quadrupole order, we use the 
SA as our first step to study the importance of fluctua- 
tions. Since derivation of this method is given in ref. 17 
we explain here some characteristics of this method and 
details of the calculation. 

In the mean-field approximation (MFA), one solves an 
effective single site problem, where an operator for inter- 
acting neighbors is replaced by an effective field. That 
field is zero in the disordered phase. On the other hand, 
the effective field in the SA is distributed with Gaussian 
distribution around a mean value even in the disordered 
phase. This fluctuating field represents the local field 
created by the short-range order. 

In order to explain the method concisely, we repre- 
sent quadrupole-moment operators as O a symbolically. 
We use the irreducible tensor operators for O a to obtain 
scalar equations from a matrix equation. cJP For H \\ (001), 
O a refers to one of the o per ators in cq^^) and/or 
eq.(||). For ff||(lll), eq-lD) and/or {-O v , - O zx + 
20 xy )/V6, {O yz - O zx )/V2 and (O yz + O zx + O xy )/V3. 
Then, the Hamiltonian for our model can be expressed 
as 

1 " 

H = 2 ? ^ J- 3 O a ^O a {j) + E Wz(i), (3.1) 

ij a — 1 i 

where n refers to the number of fluctuating components 
and is equal to 2, 3 and 5 for the r 3 -, the IV and the 
(F 3 + T5)-type interactions respectively. 



We write the mean field as a a (i) — J2j J"j(O a (j)) and 
the fluctuating field around it as ( a (i)- The effective 
single-site Hamiltonian Tli and the partition function Z\ 
is given by 

n 

Hi = E{-a"0 Q _ ^(p a - (O a ))} + Hz, (3.2) 



z 1= n 



8dC a (C a ) 2 
f4= expHjiJJL] ^TrexphSWi], 

a=i J2nf3J a J 

(3.3) 

where the site index i is omitted. The term ( t a (O a ) in 
the Hamiltonian plays a role to enforce = 0. The 
variance J a is determined by the self-consistency condi- 
tion 

d D q 



(27r)^l-(J«->) x £' 



(3.4) 



where D is the spatial dimensions (D — 3 for the present 
case), xl is the local (strain) susceptibility calculated 
with Tii, and Jg is the Fourier transform of the in- 
teraction: jg 



The lat- 



Ej J jh exp[-iq • (Rj - R ) 
tice constant is set to unity. This equation requires 
the consistency between the local susceptibility and its 
another expression using the inverse Fourier transform: 
Xl = A f_1 ^gXq- Note that this consistency is not 
satisfied in the MFA for a finite dimensional system. 
The integral can be performed analytically for the near- 
est neighbor interaction Ji* = —2J a E;=i cos< li even in 
three dimensions.^ 

The SA still keeps the noncommuting character of 
quantiim operators. In contrast, our previous theory for 
CcBgtLj) replaces the pseudo-spin vectors composed of 
the Pauli matrices by classical vectors. Further improve- 
ment of the present theory is that the renormalization is 
performed differently for each component J a when the 
model has a low symmetry. The spherical model for clas- 
sical spins used in rcf.[l^ cannot deal with the anisotropy 
in J a even if the symmetry is lowered by a magnetic field. 
These theories coincide with each other in the limit of in- 
finite temperature. 

At high temperatures, the variance J a in Gaussian dis- 
tribution of the effective field is almost zero, and there- 
fore it has little deviation from the MFA. As the tempera- 
ture decreases, the short-range order makes J a larger. In 
other words, the distribution of fluctuations gets broader. 
Finally it reaches the critical variance J cr at the critical 
temperature, and the transition occurs; the susceptibility 
for a certain component A (or for multiple components 
simultaneously) diverges with the wave number q = Q 
at which takes the maximum Jq. Hence, we obtain 



(3.5) 



and J Acr = JA - (XL") -1 - 

As the number D of spatial dimensions goes to infin- 
ity, the critical variance J CI approaches zero. In short, 
the SA approaches the MFA. This is consistent with the 
nature of the theory; the original theory for the S A is the 
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next leading order theory with respect to the inverse of 
dimensions, while the MFA is the leading order theory. 
Small J cr leads to a large transition temperature because 
the effect of fluctuation is small. On the contrary, when 
the number n of fluctuating components increases, the 
critical temperature becomes smaller. The reason is that 
the local susceptibility xl decreases with increasing n as 
shown below. 

To perform the multiple Gaussian integral, the polar 
coordinate in the C space is convenient. The reason is 
that in zero field some integrands of angular integral are 
constants by symmetry and remain smooth even in weak 
magnetic fields. For example, in terms of the polar co- 
ordinate for the (T3 + T5)-type interaction, we obtain 



n 



1 



par , jo* 



dflo I Cr d C-exp[- j. 
z 



(3.6) 



where Cr = y J2 a P(( a ) 2 /J a an d refers to the gen- 
eralized solid angle. To perform the integral on the unit 
hypersphere in five dimensions, we introduce the vari- 
ables 9 a , Ob - ■ ■ as follows: 



Ci = 

C2 = 

Ca = 

a = 

c 5 = 



Ji//3 Cr sin 9 a sin 9b sin 9 C sin 9 a , 
J2/ f3 Cr sin 9 a sin 9b sin 9 C cos 9d , 
J3 / (3 Cr sin 9 a sin 9b cos 9 C , 
J4//3 Cr sin 9 a COS 9 b , 

J5//3 CrCOS0 a , 



(3.7) 



dtte = sin 3 9 a sin 2 9 b Sm9 c d9 a d9 b d9 c d9 dl 



(3.8) 



with 9 a , 9 b , C e [0, 7r] and 9 d € [0, 2tt). 

In zero field, our model has a high symmetry so that 
Xl and J do not depend on a. Then, J2a=iXh/ n nas 
the spherical symmetry in the C space and is equal to 
XL- We obtain 



XL 



1+ n-1 



(sinh(/3C)/(/3C))< 
(cosh(/3C)) C 



(3.9) 



l + ( n-l)£WM/M^ 
J F(n/2,l/2-f3J/2)) 

where (• • -) c stands for / °° dC • • • C"" 1 exp[-/?C 2 /(2 J)] 
and F(a,7;z) is Kummer's confluent hypergeometric 
function. In fact, this equation is common to the re- 
sults of the SA for the Heisenberg model(rt = 3), the XY 
model(n = 2) and the Ising model(n = 1) in zero field. 
Let us take such spin models to see the dependence on 
n. The first term 1 in the bracket of eq.(3.10) represents 
contribution from a fluctuating field parallel to the spin 
and the second term with factor n— 1 perpendicular to it. 
We note that the quotient F{- ■ ■) / F{- ■ ■) in the second 
term of eq.(3.10) is less than unity except for f3J = 0. 
The quotient is reduced to a simple form (1 + /3J) -1 for 
n = 3, and (3 + /3J)(3 + 6/3J + /3 2 J 2 )" 1 for n=5. For 
n = 2, it takes a rather complicated form including the 



incomplete gamma function. As n becomes larger, the 
number of perpendicular components increases and thus 
the perpendicular part in eq.(3.10) becomes dominant. 
Hence it leads to the decrease of xl- It can be shown by 
more detailed investigation that xl is a monotonically 
decreasing function of n at a fixed temperature. Physi- 
cally, the short-range order of the z component disturbs 
orderings of x and y components. This corresponds in 
the SA to the fact that the distribution of C z makes x£ 
and Xl decrease and the or derin gs of x and y components 
become less favorable. In §^2, we investigate the effect 
in more detail. This effect becomes large as the number 
of fluctuating components increases. We note that the 
Ising model with n = 1 does not have such fluctuation 
effect. 

Once a magnetic field is applied, the model has a lower 
symmetry and each J a can have different values. As 
temperature decreases, on e o f them reaches the criti- 
cal value and satisfies eq^SJ)). Hence one must solve 
the set of nonlinear equations to determine the other 
J°"s and the critical temperature. We used the Newton 
method to solve the set of equations. In solving it, wc 
took care of the fact that (x£) _1 + J a m the integral 
in eq.(3.4) is always larger than Jq. At finite magnetic 
fields, the multiple Gaussian integral could not be per- 
formed analytically, and therefore we did numerical in- 
tegration. For the 1^3- and the Ts-type interactions, the 
integral is no more than a double integral by symmetry, 
and thus executable. However, for the (r3 + Ts^type 
interaction, the integral on the unit hypersphere repre- 
sented by eq.(3.8) was too time-consuming to use the 
standard routine of numerical integration. Therefore, 
we approximated the integral by summation over ten 
points (±1, 0, 0, 0, 0), (0, ±1, 0, 0, 0), • • • , (0, 0, 0, 0, ±1) on 
the unit hypersphere in five dimensions. This method 
is exact in zero field. We have checked that the same 
approximation for the I^- and the Ts-type interactions 
gives only a few percent of deviation in the transition 
temperature from the results with use of the double in- 
tegral. 

For calculation of the specific heat, we use the formula 
C = dU /dT where the internal energy U = (H) is dif- 
ferentiated numerically. In this paper, we restrict the 
calculation only to the disordered phase. Then, expec- 
tation values of quadrupole moments (O a (i)) and dipole 
moments (M(i)) do not depend on the site index i. The 
internal energy is written within the SA as 



U 



rp n 



A- lJ2 J ?A° a ) 2 - NH -( M ) 



(3.11) 



T 



Q=l 



; J^ (O a ) 2 + H ■ (M) 



(3.12) 

where \q = Xl{1 — {Jq ~ ^ q )Xl} _1 and the site index 
is omitted. We used here the relation A -1 J2 q JqXq = 



J°Xlj which is derived from eq.(3.4). All of the quan- 
tities in eq.(3.12) can be calculated with the effective 
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single-site Hamiltonian Hi. 

§4. Numerical Results 

4-1 Quadrupolar ordering temperature 

In contrast to the MFA, the present effective medium 
theory with the SA gives the result that Tq(H) increases 
with H . Figure 1 shows the phase diagrams obtained by 
the present theory and those by the MFA for comparison. 
The low temperature phase is the antiferro-quadrupolar 
phase. The magnetic field is oriented along either (001) 
or (111). As magnetic field increases, numerical integra- 
tion of Gaussian distribution becomes difficult. Hence 
the phase diagrams are plotted up to the maximum field 
we can calculate. 

We take the interaction constant J' as the unit of en- 
ergy, where J' refers to either Jf 3 or Jf 5 in eq.(2.10). 
Accordingly J'/fce is taken to be the unit of tempera- 
ture and J'/mb the unit of magnetic field. In this unit, 
T Q (H = 0) in the MFA is 6 (= z n ) irrespective of the 
number of competing order parameters. With inclusion 
of fluctuations we have the following relation: 



At finite H 1 Tq(H) shows some increase except for 
the T3-type with ff||(lll). The rate of increase R = 
{max[T Q (ff)] - Tq(0)}/T q (0) satisfies i? r3 < R r5 < 
i? r3+rs . That is to say, the rate becomes larger as the 
number of components of interaction increases. 

The types of order parameters on the phase boundary 
are summarized in Table I. A cusp in the phase diagram 
means the change from (O2 ) to (O2 2 ) with increasing 
H. The magnetic field H c at the cusp is also shown 
in Table I. The phase boundary for the (r 3 + r 5 )-type 
interaction in weak fields has very small dependence on 
the orientation of magnetic field. Although Tq(H) in 
weak fields seems to have almost linear field-dependence 
in Figure 1, it is in fact perpendicular to the T-axis at 
H — 0. Since we are interested mainly in the vicinity of 
the phase boundary between the phases I and II, study 
of the order parameters at T < Tq(H) is left for future 
study. 

4-2 Specific heat and entropy 

Figure 2 shows the specific heat for the case of the 
Ts-type interaction with ff||(lll) above the transition 
temperature Tq(H). We notice that the MFA gives zero 
specific heat in zero field above Tq (H) . It is known in the 
SA for the Ising model that the specific heat has a cusp 
at the transition without a jumpJla) We expect a similar 
behavior in our model although we have not calculated 
the specific heat below the transition temperature. At 
weak magnetic field the specific heat decreases monotoni- 
cally as temperature increases above Tq (H) . It is evident 
that the specific heat near Tq(H) becomes sharper and 
larger with increasing magnetic field. With further in- 
crease of the field beyond [XbH ~ 5 J', the peak becomes 
smaller and the Schottky type anomaly appears. This is 
due to the Zeeman splitting of the r 8 levels. We discuss 
implication of the result to experiment in §[| 

We calculate the entropy relative to the value at the 



transition temperature by numerical integration of the 
specific heat. Figure 3 shows the result. In the high 
temperature limit the entropy should tend to In 4. Hence 
we can roughly estimate the value at Tq(H = 0) since 
it almost saturates in the high temperature end of the 
calculated range. It is remarkable that the entropy at the 
transition temperature is much less than the full value 
In 4 even in the absence of Zeeman splitting. We will 
report the details of thermodynamic quantities including 
magnetization for various direction of H and different 
types of interactions in a separate paper. 

§5. Mechanism of Increased Transition Temper- 
ature by Magnetic Field 

5.1 Decreasing number of fluctuating components 
Our results show pronounced difference of the phase 

diagram and entropy between the MFA and the SA. The 
simplest model we take proves to be a convenient model 
to investigate the effect of the fluctuations. In this sec- 
tion we discuss how a magnetic field suppresses the fluc- 
tuation. 

When different parameters are competing for the sta- 
bility, the transition temperature should be lowered in 
general. If there are many competing order parame- 
ters, different types of short-range order should occur 
and disturb each other. Therefore the transition tem- 
perature should become lower than that of a system in 
which only one of them is favored. For our model, the 
r 3 -type interaction has two competing components, the 
T5-type has three and the (r 3 +T5)-type has five. Then, 
Tq(H = 0) in the SA decreases with increasing number 
of components. The reason is that a model with many 
components has many types of competing orderings. On 
the other hand, the MFA cannot deal with this compe- 
tition and gives the same Tq(H = 0) for all types of 
interaction. 

This type of competition should be suppressed with in- 
creasing H in our model. In other words, the ordering of 
the one (or two) of several components becomes more fa- 
vorable than the others, and the number of competing or- 
derings decreases. Consequently Tq(H) should increase. 
This effect becomes larger as the number of components 
increases because the situation changes more drastically 
when one of a larger number of components becomes fa- 
vorable. Therefore, this mechanism naturally explains 
the result that the rate R of increase in Tq(H) satisfies 
i? r3 < R r5 < i? r3+rs . On the contrary, we obtained 
monotonically decreasing Tq(H) for the r 3 -type inter- 
action with In this case, such suppression of 
fluctuations does not occur because the ordering of O2 2 
is as favorable as that of O2 owing to the trigonal sym- 
metry. In order to confirm this mechanism, we calculated 
Tq(H) in another model which has only one component 
like the Ising model. This model does not have such 
mechanism of suppression. We indeed find that Tq(H) 
does not increase with H . 

5. 2 Effects of fluctuating magnetic field on spin order- 
ing 

The SA applied to the Heisenberg model and the XY 
model does not give increasing transition temperature 
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with increasing magnetic field. One difference between 
the Hcisenberg model and the f^-type interaction model 
is that while only one component among three compo- 
nents becomes favorable at finite magnetic fields in the 
latter, two components perpendicular to the magnetic 
field become favorable in the former. As a result, the sup- 
pression in the Heisenberg model is rather small. How- 
ever, this remark cannot be applied in comparing the 
XY model and the I^-type interaction. Another dif- 
ference, which is probably more important, is that a 
magnetic field is not conjugate to quadrupole moments 
and is coupled with them in a complex way. To ana- 
lyze this, it is convenient to use the pseudo-spin repre- 
sentation of the Zeeman term with ff||(001), namely, 
{a z + (4/7)t z ct z }[1bH . The quadrupole moment t z 
(O2 ) couples to a kind of fluctuating field a z which has 
the Ising-type distribution. As shown below, the fluctu- 
ating field a z has a character of suppressing fluctuation of 
quadrupole moments. However, the Zeeman term plays 
a dual role of favoring a particular type among compet- 
ing fluctuating fields, and also destroying the quadrupo- 
lar ordering as in antifcrromagnets. In the following, we 
first take a spin model to clarify effects of a fluctuating 
magnetic field. Then, we investigate how a static mag- 
netic field affects the fluctuation in the I^-type interac- 
tion model. It is shown that the dominant effect of the 
Zeeman term at weak field is suppressing quadrupolar 
fluctuations. 

First of all, we take a single classical spin coupled with 
a Gaussian fluctuating field. The Gaussian identity, 

f d f expH3^£wS'] = exp[/?— (S z )% 
^ y/2nfjV z 1 2V Z A y 1 2 J ' 

(5.1) 

shows that a Gaussian fluctuating field ip z coupled to a 
classical variable S z is equivalent to a field V z jl coupled 
to (S z ) 2 . That is to say, ((S z ) 2 ) increases with increas- 
ing V z and the relation -2dVL/dV z = {{S z f) is satis- 
fied where Q is the thermodynamic potential. Therefore, 
the fluctuating field increases the susceptibility \ z — 
P(((S Z ) 2 ) — (S z ) 2 ). On the contrary, x x and \ v are de- 
creased because of the identity {(S x ) 2 + (S y ) 2 + {S z ) 2 } = 
S(S +1). Note that the increase of x z , and the de- 
crease of \ x an d X v are °f 0{V Z ) because V z is a field 
conjugate to (S z ) 2 . Even for a quantum spin, the discus- 
sion above is still valid provided that S z is replaced by 
Jq dr5 z (r). Namely, the susceptibility \ z is given 

by 



as 



drS z (T) 



dr(S z (r)S z (0)). 



(5.2) 



(5.3) 



A difference from a classical spin is that the suscep- 
tibilities are not directly determined by the identity 
({S x ) 2 + {Sv) 2 + (S z ) 2 ) = S(S + 1). For the case of 
S = 1/2, the susceptibilities in the presence of the Gaus- 
sian fluctuating field with variance V z / j3 are calculated 



P F(l/2,'3/2;PV z /2) 
X 4' X 4F(l/2,l/2;/3y z /2)' K ' 

X x has the same form as the quotient in the second term 
of eq.(3.10) and is a monotonically decreasing function 
of V z . Therefore, the fluctuating field decreases x x an d 
X v , but does not affect x z ■ The reason for the latter is 
that the Hamiltonian and S z are commutative, and that 

one has j/?" 1 f£ dTS z (r)} = (S z ) 2 = 1/4 independent 

of V z . For S > 1/2, (S z ) 2 is not a constant and thus 
X z is increased by the fluctuating field. Note that in any 
case the change of susceptibilities is of first order in V z 
by the same reason as in the classical system. 

Next, we treat the lattice systems. Let us take the 
following model as an example which have both the 
Heisenberg-type exchange interaction and Gaussian fluc- 
tuating fields: 



H 



Jii S 'i ■ Si 



2 

i'j 



(5.5) 



For classical spin systems and quantum systems with S > 
1/2, the local susceptibility xt at each site is increased 
by the Gaussian fluctuating field as well as that of the 
single spin. Consequently, the transition temperature to 
the ordered state should increase. In the mean-field ap- 
proximation, the susceptibility xq = Xl/(1 — JqXl) for 
a certain wave number Q diverges at the phase transi- 
tion. Hence the transition temperature increases if xt i s 
increased by the fluctuating field iff. Even in quantum 
spin systems with S = 1/2, the transition temperature 
should still increase by the fluctuation effect. As V z 
increases, ordering of S x or S v should become less fa- 
vorable. Therefore, in the limit of large V z , only the S z 
ordering should be possible. This means that the fluc- 
tuation among components is suppressed. Certainly, the 
transition temperature increases with V z in the SA to 
this model. Further, the local susceptibility xt increases 
here. The reason is that V z makes xt an d Xl decrease. 
Then, J x and J y decrease and xt increases. On the 
other hand, when we solve the model by the mean-field 
approximation, the transition temperature T c for the S x 
and S y orderings decreases with increasing V z , while T c 
for the S z ordering does not change. Namely, in the spin 
model defined by eq.( |5.5|) with S — 1/2 the transition 
temperature does increase with field in contrast to the 
result of the mean field approximation. 

This fluctuating field has similar effect even if ip z is 
not a Gaussian-distributed variable but the Ising-type 
variable taking values of ip z = ±^fV z / (3. This field re- 
sembles ha z in the Zeeman term with if||(001) repre- 
sented by the pseudo-spins for the quadrupolar model. 
We discuss it in the next subsection. 

5.3 Two competing effects of magnetic field on the 
quadrupolar ordering 
In the pseudo-spin representation, the interaction 
Hamiltonian for the T3-type interaction is the same as 
that of the XY model. The only difference between them 
is the Zeeman term. Nevertheless, their phase diagrams 



Fluctuation Effects on the Quadrupolar Ordering in Magnetic Field 



V 



by the SA with J?||(001) show a remarkable difference; 
while the transition temperature increases with H in the 
I^-type interaction, it does not in the XY model. In 
this subsection we discuss the reason for the difference 
in weak fields. 

The Hamiltonian of the r 3 -type interaction with term — (4/7)(3h 2 Tf to the strain susceptibilities should 
fT||(001) is written as 



the Zeeman term in the XY model and thus it works to 
induce (rf) uniformly and to destroy the antiferro-type 
ordering. Since the lowest order contribution of a mag- 
netic field to the susceptibilities in the XY model is the 
square of the field, the lowest-order contribution of the 



n = \ E J vW + t * t D - E p* Ha *i ( 1 



4 . 



H 



int 



ha? 1 



4 . 



(5.6) 
(5.7) 



where h refers to (IbH. In this expression, the er-spin 
emerges only as a z . Hence, one can treat it as a classical 
variable taking values ±1. To study the structure of the 
partition function, we introduce path-integral represen- 
tation for the r-spin. Since the explicit expression for 
the Berry phase term is not necessary here, we write it 
symbolically as £b- 

Drexp -/ dr{£ B (r) +H(t)} , (5.8) 



E 







be the square of -(4/7)/3/i 2 , namely 0[(/3/i) 4 ]. 

On the other hand, the Gaussian fluctuating field 
<p* with variance (4/i/7) 2 affect the susceptibilities by 
0[((3h) 2 ] as described in the previous subsection. There- 
fore, it should make a dominant contribution at weak 
fields. This term makes the ordering of r z favorable. Af- 
ter all a weak magnetic field suppresses the fluctuation 
and increases the transition temperature. 

In the previous paper of Uimin and the present au- 
thorsJifP the pseudo-spins for the quadrupole moments 
are approximated by classical spins. In the previous sub- 
section, we have shown that the fluctuating field ip z has 
smaller influence on \ z m the quantum spin with S = 1/2 
than that in the classical spin. Hence the approximation 
in ref.[t5] enhances the effect of the fluctuating field iff in 
eq.(5.12). The classical treatment should be the reason 
why the result of the IVtype interaction with H ||(001) 
shows much larger increase of Tq(H) than that of the 
present paper. 



After summation over a z , we neglect terms of higher 
order in (3h to investigate the behavior in weak fields. 
We use the approximation cosh a; ~ e x / 2 for x <C 1. 
Here we have 

> h I ( l + y7f(r)^dr. 



Then we use the identity eq.(5.1) to introduce the fluctuating field iff. Finally we turn to the operator representation 
with Tr instead of J T>t. The calculation explained above goes as 



T>t exp 



Vrft 



{C B {T)+H int (T)}dT 




(5.9) 



d<pf 



-oo V2nh 2 



exp 



{C B (T)+H in t(T)}dT 



2 -n 
2 *n 



V2^P 
7d^f 



E 



(vf) 2 



2h 2 



4:V2nh 2 



exp 



exp 



E 21, 1 

Y. 1 lvf> 



Tr exp 



& f (l + ^7f(r))dr 



(5.10) 



(5.11) 



2{Ah/7f 



-B 2 h 2 



x Trexp 



■0 {H, 



E 



(5.12) 



In the expression (5.11), the only difference from the 
original expression (5.6) is that the Ising-type distribu- 
tion ha* is replaced by ip*, which obeys the Gaussian 
distribution with the variance h 2 . Then, we replaced (pf 
by (7/4)</?f + (3h 2 in deriving eq.(5.12). 

The linear term with respect to h is absent in eq.(5.12) 
and t z couples to h 2 because of the time reversal sym- 
metry. The term -(4/7)/3/i 2 rf is of the same form as 
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§6. Discussion and Concluding Remarks 

In this paper we have shown for quadrupolar interac- 
tion systems that a magnetic field has two contrasting 
effects. One effect is inducing quadrupole moments uni- 
formly and destroying the antiferro-quadrupolar order- 
ing. The induced quadrupole moments are even func- 
tions of magnetic field owing to the time reversal sym- 
metry, and the change of Tq(0) is very small (oc H A ) at 
weak fields. The other effect is suppressing the compe- 
tition among different components of order parameters. 
Since the resultant change is of 0(H 2 ), the latter effect 
dominates at weak fields. Hence the transition tempera- 
ture increases as H increases. 

For comparison with experimental results in CeBg, a 
crude estimate of J' is obtained by comparing the theo- 
retical value of ?q(0) with the experimental one. With 
the value of J' so determined, the unit for magnetic field 
is also fixed. Let us assume for the moment the T§ or- 
der. Since the numerical value of Tq(0) in units of J' /fcs 
is the same order as the experimental value in units of 
Kelvin, one may roughly estimate J'/fce ~ IK, and the 
unit for magnetic field as 1.5 T (tesla). Then the in- 
crease of Tq(H) in Figure 1 turns out not so strong as 
the experimentally observed one, and the re-entrant field 
is much smaller than the actual one which is larger than 
the maximum of available magnetic field. 

On the other hand, results of the specific heat at weak 
fields agree qualitatively with experimental ones both for 
the magnitude and the temperature dependence. We 
emphasize that the lowering of the cubic symmetry is 
not necessary to explain the entropy smaller than In 4 
at the transition temperature, and the large tail of the 
specific heat experimentally observed. The temperature 
range of a weak Schottky type anomaly found in Figure 
2 for if ~ 12 T and larger overlaps actually with the 
experimental Tq(H) up to 15 T in CeBg.EP Hence the 
anomaly is not seen experimentally. However in the alloy 
system Geo.5Lao.5B6, such Schottky anomaly is indeed 
observed^ due to the lowered Tq(H). 

For serious comparison with experimental results, we 
have to consider two directions of development on the 
basis of the present study. The first is to choose a realis- 
tic model to describe the interaction among quadrupole 
moments at different sites within the SApConcerning the 
intersite interaction, the rplpjaLtbe rangeH^ and types in- 
cluding higher multipolcslla'LJ'ES should be reexamined 
with inclusion of fluctuations. With the intersite octupo- 
lar interaction, the effect of magnetic field on the shift of 
Tq(H) becomes larger and should compare more favor- 
ably with experiment. Although we have neglected the 
higher crystal field levels other than the ground level, the 
crystal field splitting is so small in TmTe that the Tq and 
the doublets seem to be within a few meV above the 
ground r 8 level. Hence, it is interesting to investigate 
how the crystal field splitting affects the phase diagram 
of the quadrupolar ordering. 

Another direction of development is including dynam- 
ical and quantum effects of fluctuations. The most im- 
portant may be the Kondo effect in the case of CeBg but 
not in the case of insulating systems such as TmTe. In 



the latter case the quantum fluctuation of multipolc mo- 
ments should persist even without the Kondo effect S3 
In order to deal with dynamical effects in the effective 
medium theory, numerical methods such as the resol- 
vent method, the quantum Monte Carlo or the numeri- 
cal renormalization group seem promising. We hope to 
report on these development in the near future. 
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Table I. Order parameters on the phase boundary. Here H c refers to the critical magnetic field at which the type of order parameters 
changes. 
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Fig. 1. The phase diagrams obtained by the SA(solid lines) and the MFA(brokcn lines). 




Fig. 2. The specific heat above the transition temperature obtained by the present theory for the Ts-type interaction with if||(lll). 
Each transition temperature is indicated by an arrow. The unit of the ordinate corresponds to 8.31J/(mol- K). 




